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In order to treat tensor force explicitly, we propose a microscopic model for nuclear 
structure based on antisymmetrized molecular dynamics (AMD). As a result of the present 
study, it is found that some extentions of the AMD method are effective to incorporate the 
tensor correlation into wave functions. Calculating the deuteron, triton and 4 He with the 
extended version of AMD, we obtained solutions that have large contribution of the tensor 
force. By analyzing the wave funtion of 4 He, it is found that its single-particle orbits are 
composed of normal-size Os orbits and shrunk p orbits so as to gain the tensor force. 

§1. Introduction 



Tensor force plays important roles concerning various properties of nuclei. Re- 
cent studies with oh initio calculations show that the contribution of the tensor 
force to binding energies is large in A ~ 10 light nuclei, 1 ) as well as in few-nucleon 
systems. 2 ^ Besides the binding energies, effects of the tensor force on single-particle 
energies and nuclear structure have been suggested for a long time. 3 )' 4 ) Recently, 
several studies showed its influences on LS splitting, 5 ) GT transitions and magnetic 
moments. 6 ) These facts imply that direct effects of the tensor force may be revealed 



in some properties of nuclear structure. 



The tensor force causes a characteristic correlation between two nucleons, which 
is called the tensor correlation. In the deuteron, the tensor force gives a large attrac- 
tion due to the coupling of 3 Si with 3 Di. Also, in 4 He, the [JLS = 022] component 
is mixed to the dominant (0s) 4 component due to the tensor force. Since the tensor 
correlation between two nucleons is varied, depending on the environment surround- 
ing them, the contribution of the tensor force is sensitive to the nuclear density and 
structure. According to analyses with the g matrix, as the density increases in nu- 
clear matter, the contribution of the tensor force is suppressed (tensor suppression). 
This is a major cause of the nuclear saturation property. 7 ) Reversely, as the density 
decreases, the tensor suppression becomes small. Therefore, the tensor force is con- 
sidered to be important in light nuclei because they have a relatively large surface 
region, being low density. In Ref. [4], it is emphasized that the renormalization 
of the tensor force to an effective 3 E central force is very sensitive to the starting 
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energy as well as the density. Because the starting energy is strongly dependent on 
the structure of the total system, the tensor force gives a structure dependence to 
the effective interaction. Thus, the tensor force might be closely related to a variety 
of structures realized in light nuclei such as a clustering structure. Moreover, it is 
natural to expect that the tensor force may provide clarity to the problems of un- 
stable nuclei, which usual approaches with effective nuclear forces sometimes fail to 
describe. 

There are many kinds of model calculations that are useful for systematic stud- 
ies of various nuclei in wide mass-number regions. They, however, have avoided any 
explicit treatment of the tensor force because of its difficulty. In such model cal- 
culations, one uses effective interactions that are either derived from realistic forces 
theoretically by renormalizing the tensor force into central and LS forces, or phe- 
nomenologically determined. In usual cases, the density dependence arising from 
the tensor force is simulated by phenomenological density-dependent or three-body 
forces. Such model calculations with effective interactions have succeeded to repro- 
duce and explain many physical quantities of various nuclei. However, there remain 
some phenomena that those model calculations failed to reproduce without using 
system-dependent effective forces. The parity inversion problem of 11 Be and the 
binding mechanism of 11 Li are typical examples. Therefore, it is natural to expect 
that the effect of the tensor force might be a key to solve these problems. 8 ) More- 
over, another drawback of the usual model calculations is that one can not discuss 
the effect of the tensor force on physical quantities in an explicit way, because the 
obtained wave functions do not contain a tensor correlation. 

Our aim is to treat the tensor force as directly as possible in the model calcula- 
tion. We would like to systematically investigate its contribution to the binding of 
various nuclei and the effects of the tensor correlation on the nuclear properties. We 
mention here a repulsive core of the nuclear force, which is another characteristic of 
a realistic nuclear force, and is not explicitly contained in the usual effective nuclear 
force as well as the tensor force. Recent developments in treating the repulsive core 
and the tensor force in a model calculation can be seen in the unitary correlator 
method, 9 -* the ^-matrix method 10 ) and the Vi ow _k method. 11 ) Since the repulsive 
core should be less state-dependent than the tensor force, we consider that it can 
be smoothed out for the present purpose. We adopt an effective force where the re- 
pulsive core in the central force is smoothed by the ^-matrix method in the present 
study. 

To accomplish our purpose, we need a wave function that can contain the tensor 
correlation. We start from a mean-field approach, considering its applicability to 
wide mass-number regions. As already known, ordinary mean-field wave functions 
are insufficient to treat the tensor force. From the operator form of the tensor force, 
one can easily expect that two extensions are needed. The first one is i) a parity- 
violating mean field and a strong correlation between spin and space. The other point 
is ii) flexibility of the isospin wave function. The tensor force contains a parity-odd 
operator as a dominant term in a mean-field picture, and it also includes a charge- 
exchange term. Therefore, it is difficult to treat the tensor force within the usual 
mean-field approaches. In fact, the expectation value of the tensor operator vanishes 
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in a model space of a simple mean-field. The importance of the parity-violating 
mean field and charge mixing are already pointed out in the Ref. [12]. 

Several groups are currently challenging this theme. Sugimoto et al. are investi- 
gating the contribution of the tensor force in various nuclei, based on the relativistic 6 ) 
and non-relativistic 12 ) mean-field approaches. Myo et al. are investigating the effect 
of a tensor correlation in light unstable nuclei with an extended shell model. 5 ) We, 
here, adopt an antisymmetrized molecular dynamics (AMD) method 13 -* as a starting 
point of our model. Due to the simplicity of its framework, AMD can be easily 
improved, while corresponding to our requirements. 14 )~ 17 ) Since its wave function 
is simply composed of Gaussian wave packets, we can easily perform various pro- 
jections, such as parity, angular momentum and charge number projections. The 
parity- violating mean field has already been realized by the parity projection in 
AMD. Other advantages of the AMD are that it does not rely on any model assump- 
tions on nuclear structure, such as the existence of clusters and nuclear deformation, 
and that it is applicable to general nuclei over a wide mass-number region. It has 
succeeded to describe various structures and cluster aspects of light and middle nu- 
clei. 13 ^ 17 ) We, therefore, expect that a model based on AMD can directly reveal 
a relation between the tensor force and various structures of light nuclei, including 
clustering structure, as investigated in Ref. [4]. On the other hand, in Refs., 5 ^ 6 -*' 12 -* 
the model wave functions are based on mean-field approaches, and the single-particle 
orbits in 4 He are assumed to be linear combination of s and p-orbits. Based on the 
AMD method, we will see how the tensor correlation is realized in the wave functions 
without relying on such a model assumption. 

In the present paper, we employ effective interactions that involve the tensor 
term explicitly, and extend the framework of AMD so that it can gain the contri- 
bution of the tensor force. We adopt an effective force that is derived from the 
Tamagaki potential 18 ^ by smoothing the repulsive core of the central term, and also 
test another effective force: the Furutani- Tamagaki potential. 19 ) We calculate the 
deuteron, triton and 4 He with an extended version of AMD. We show the effects 
of i) a parity-violating mean field and a strong correlation between spin and space 
and ii) the flexibility of the isospin wave function in practical calculations. More- 
over, we propose additional extensions of the AMD to gain the tensor force more 
effectively; iii) linear combination of the wave packets that have different width pa- 
rameters and iv) charge-, parity- and angular-momentum-projections. Finally, we 
analyze the properties of the wave functions of the deuteron and 4 He obtained by 
the extended AMD, while focusing on the tensor correlation. 

This paper is organized as follows. In section 2, we explain the extensions of 
the AMD framework. We investigate the effect of these extensions in the deuteron, 
triton and 4 He in section 3. Analysis of the obtained wave function of 4 He from 
the viewpoint of single particle levels is shown in section 4. In section 5, we give a 
summary and some discussions on further developments. 



4 A. Dote, Y. Kanada-En'yo, H. Horiuchi, Y. Akaishi and K. Ikeda 



§2. Formalism 

2.1. Hamiltonian and tensor force 

The Hamiltonian used in this paper is 

A A 

H = E + E ( u £ + 4 + ^ om6 ) - 7cm. (2-1) 

i=l i<cj=l 

nr\2 

Here, U is a kinetic energy term, (p i , momentum operator; m^, nucleon mass); 
vfj, vfj and v 9? ulomh indicate the central force, tensor force and Coulomb force, re- 
spectively. The center-of-mass motion energy, Tqm = 2Am N {^2t=i Pi) > * s removed. 
We employ two kinds of effective interactions, the Furutani-Tamagaki potential (F- 
T potential) 19 ** and the Akaishi(T) potential (AK(T) potential), as a set of vfj and 
vfj. Both potentials are described by a superposition of several range Gaussian wave 
packets: 

N 

v?2 = E { C l Ep CE) + Cl E P?E) + C^PCO) + C 3 n °P( 3 0)} exp [-(r/6 n ) 2 ^2,-2) 

n=l 
N 

vj 2 = E Tn E P( 3 E) r 2 S 12 exp [-{r/b n f} , (2-3) 

n=l 

where r = \n - r 2 \. N = 3 in the F-T potential and N = 10 in the AK(T) 
potential. P(X), X = 1 E, 3 E, 1 0, 3 0, is a projection operator and 5i2 is a tensor 
operator, S12 = 3(<Ti • t)(<T2 • r)/r 2 — (<ri • (T2), where r = r\ — r 2 - The radial 
part of the tensor force in both potentials is represented by a linear combination of 
r 2 x Gaussians. The interaction parameters, {C* ,T 3E , b n }, in the F-T potential are 
given in Ref. [19], and those in the AK(T) potential are given in Appendix. The 
Coulomb force, v f. oulomh ^ \ s expressed by the superposition of seven range Gaussian 
wave packets. 

The F-T potential was used in a cluster model study of 3 He+p scattering, and 
succeeded in reproducing scattering data, including polarization quantities and the 
spectra of 4 Li. 19 ) It can also reproduce the binding energy of 4 He with a harmonic- 
oscillator (0s) 4 wave function. The AK(T) potential is derived from the Tamagaki 
potential (OPEG), 18 ) which is a realistic nucleon-nucleon potential. The central part 
of the AK(T) potential is obtained by smoothing the repulsive core of the central 
part of OPEG with the method of ^-matrix theory, whereas the tensor part remains 
to be the same as the original one of OPEG. In other words, the tensor force is 
not renormalized into the central one when constructing the AK(T) potential. The 
details of its derivation are explained in Ref. [10].*) The 1 E and 3 E central forces 

*-* In Ref. [10], the AV8' potential 21 ^ is used and several values of the cut-off momentum, Kq \ 
for the tensor force are investigated. The AK(T) potential corresponds to the case that the Tamagaki 
potential is used instead of the AV8' potential and fcg 5 ' = 00 fm -1 . 
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Fig. 1. (a) Furutani-Tamagaki potential and (b) Akaishi(T) potential. The 1 E and 3 E central 
potential and the 3 E tensor potential are shown as 'Central IE', 'Central 3E' and 'Tensor 3E' 
in each panel, respectively. Both potentials are multiplied by r 2 . 
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(b) Akaishi(T) potential 



and the S E tensor force in both potentials are depicted in Figure 1. In the F-T 
potential, the 3 E central force is dominant, while the tensor force is smaller. On the 
other hand, the 3 E central force in the AK(T) potential is smaller than that in the 
F-T potential. The tensor force has almost the same magnitude as the 3 E central 
force in the AK(T) potential. 

2.2. Extended framework of AMD 

The details of the usual AMD framework are given in Refs. [13] and [20]. In order 
to incorporate the effect of the tensor force efficiently, we improve the single-nucleon 
wave functions, \ipi), by superposing Gaussian wave packets, as follows: 



-v n r 



N /9 v 3/4 
a=l ^ ' 

K)=Q + ^)|T) + Q-^)II), 

i7i> = Q+7i) \p)-(l- j 



(2-4) 

(2-5) 
(2-6) 



where \(3 l a ) is an intrinsic-spin wave function and [7^) is an isospin wave function. 
Coefficients C x a , (3^, 7^ and the Gaussian center Z l a are complex variational param- 
eters. This single-nucleon wave function, \<pi), is more sophisticated than that in 
our previous study, 15 -* where only the spatial part of a single nucleon wave func- 
tion was expressed by the superposition of Gaussian wave packets, regarding several 
points. In the present wave function, not only the centers {Z l a }, but also the width 
parameters {u a } of the wave packets are different from each other. Moreover, the 
variational intrinsic-spin wave function, \(3'^), is independently given for each wave 
packet. Since the present single-nucleon wave function, \<pi), is possible to describe 
the strong correlation between the spatial part and the spin part, it can effectively 
take the contribution from the term <r • r in the tensor force, as mentioned in the 



6 A. Dote, Y. Kanada-En'yo, H. Horiuchi, Y. Akaishi and K. Ikeda 



previous section. In the isospin part, we also use the variational isospin wave func- 
tion, |7„), as shown in Eq. (2-6), so that it can contain both the proton and neutron 
components, though it is fixed to be a proton or neutron state in the usual AMD. 
Due to the flexibility of the isospin wave function, the present wave function can take 
effect of the term t\ • T2, especially t\ + T2- and ti_T2+, in the tensor force. This 
corresponds to point ii) mentioned in the previous section. 

In our treatment, the total wave function is constructed from an intrinsic wave 
function, |#), with various projections, where |#) is a Slater determinant of {|</>i}}: 
|#) = det[|c^j)]. We performed parity projection (P-p) and angular momentum pro- 
jection (Pj) in the same way as for the usual AMD calculation for a nuclear structure 
study. In addition, we need charge-projection (Pr z ) in the present framework, since 
|^>) is a charge-mixed state, as shown in Eq. (2-6). The charge projection is done by 
rotating around the z-axis in isospin space, 



Thus, the total wave function in the present study is the parity-, angular-momentum- 
and charge-number-eigen wave function, 



Our procedure of energy variation and projections is as follows. We employ 
\P Tz P-p<p) as a trial wave function in the energy variation. The complex variational 
parameters, {X l a } = {C % a , Z l a , (3 l a . 7^}, are determined so as to minimize the energy, 



by the "frictional cooling method with constraints," which is one of the energy- 
variation methods. After the energy- variation, \Pr z Pp <&) is projected onto the eigen- 
state of the total angular momentum; \PjPr z Pp<I>). We evaluate various quantities 
with the wave function \PjPr z Pp<P). 

We here comment on the practical scheme of the frictional cooling method with 
constraints. As mentioned in Ref. [22], we need a constraint condition to fix the 
position and momentum of the center of mass ((Rg) an d (Pg)) to the origin, when a 
single nucleon wave function is represented by the superposition of several Gaussian 
wave packets. In the present study, such several constraint conditions are satisfied by 
adding harmonic-oscillator- type potentials to the energy, TC, in the frictional cooling 
equation: 




(2-7) 



\PjP Tz Pr$). 



(2-8) 



n = 



(P Tz P v $\H\P Tz P v $) 

{p Tz p v $\p Tz p r $) 



(2-9) 



X t a = (\ + i l j) 



and C.C., 



(2-10) 



M 



n' = n + c 1 {(k g ) 2 + <p G > 2 } + Y,c a (w a - w° a f 



(2-11) 



a=2 



where W a and W® indicate an expectation value of some operator and a constant 
value, respectively. By adopting adequately large positive constant values of C\ 
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and C a in this equation, we can obtain an optimum solution that satisfies several 
constraint conditions {W a = W®} after frictional cooling. This method works much 
better than the Lagrange-multiplier method. 15 )' 22 ) 

§3. Practical results for deuteron and 4 He 

As explained in §2, we drastically extend the AMD wave function for treating 
the tensor force. In this section, we show the effect of each extension on the deuteron 
and 4 He. 

We consider four extensions: i) a parity- violating mean field and a strong corre- 
lation between spin and space, ii) flexibility on the isospin wave function, iii) different 
width wave packets and iv) charge-, parity- and angular momentum projections, i) 
and ii) are attributed to the operator form of the tensor force, (r • r) (<Ji ■ r)(<T2 • r). 
Extension i) is effective to the operator <Tj • r, because this operator changes the 
parity of a single-particle state and leads to a strong correlation between spin and 
space. In the AMD method, the former part of i) has already been realized because 
the parity projection is usually performed. The latter part of i) is taken into account 
in the present framework by the superposition of Gaussian wave packets that have 
different spin wave functions. The other operator, t\ ■ T2, causes charge- mixing in 
the single-particle wave function, which is treated by extension ii). We first examine 
the influence of the tensor force on the properties of the wave function with these ex- 
tensions. In other words, we look into the perturbative change from the usual model 
calculations obtained by an effective central-type interaction. Therefore, we use the 
F-T potential, where the tensor term is added to a central interaction, because such 
properties as the binding energy and radius of 4 He are reproduced with the central 
term in the F-T potential by simple AMD calculations without any extensions. Next, 
we adopt the AK(T) potential, which is derived by smoothing out the hard core of 
the central term in the ^-matrix theory based on a realistic nucleon-nucleon poten- 
tial. In practical calculations with the AK(T) potential, we show that extensions iii) 
and iv) are important to describe the properties of 4 He. 

In §§3.1 and §§3.2, we investigate the effects of extensions i) and ii) with the F-T 
potential, in a simplified model space where the width parameters of all wave packets 
are set to a common value, namely v a = v. The effect of extension iii) by using 
different width parameters is investigated in §§3.3, with the AK(T) potential. In all 
calculations, the deuteron and 4 He are projected onto J T = 1 + and + , respectively. 

3.1. Superposition of Gaussian wave packets with different spins and parity projec- 
tion 

As mentioned before, in order to treat point i) a parity-violating mean field 
and strong spin-space correlation in the AMD framework, the parity projection and 
the superposition of wave packets that have independent spins are needed. We 
investigated their effects by calculating the deuteron with the F-T potential. Fig. 
2(a) shows the binding energy and the expectation value of the tensor force (Vt) 
for various numbers "iV" of wave packets superposed in \ipi). Here, we have fixed 
each isospin wave function to be a proton or neutron; \ j l a ) = \p) or \n) in Eq. (2-6). 



8 A. Dote, Y. Kanada-En'yo, H. Horiuchi, Y. Akaishi and K. Ikeda 



Fig. 2. (a) Binding energy and expectation value of the tensor force (Vt) in the deuteron as a 
function of the number (N) of wave packets superposed to compose a single nucleon wave 
function, as described in Eq. (2-4). (b) Expectation value of squared orbital angular momentum 
in deuteron for various N. The solid (dashed) lines indicate the results with independent 
(common) intrinsic spins for the wave packets in each single nucleon wave function. Here, the 
F-T potential is used for both (a) and (b). The details are explained in the text. 




(a) (b) 



"common spin" (dashed line) indicates that the common spin wave function is used 
for all wave packets in each \<pi), namely \(5V) = \j3 l ) in Eq. (2-4), whereas "different 
spin" (solid line) indicates that each wave packet independently has a different spin 
wave function, as shown in Eq. (2-4). 

As shown in the figure, the expectation value of the tensor force, which we de- 
note as \Vt\, is nearly equal to zero when the superposition of wave packets is not 
performed (N = 1). In the case that different spin wave functions are used for each 
wave packet, \Vr\ is increased and the binding energy is gained by the superposition 
(see the N >2 region of "different spin" in Fig. 2(a)). A drastic energy gain of the 
Vt term is seen when the number of wave packets (N) is changed from 1 to 2. On 
the other hand, when we use a common spin wave function for each wave packet, 
\Vt\ remains to be nearly zero, even if more than two wave packets are used (N > 2 
and "common spin"). In order to see the mixing ratio of high-angular- momentum 
components, we show the expectation value of the squared orbital angular momen- 
turn, L 2 = (L ), of the total system in Fig. 2(b). Corresponding to the behavior of 
Vt, L 2 is also drastically deviated from zero in the case that a different spin wave 
function is used (solid line) . This deviation of L 2 from zero indicates a mixing of the 
D-state, which is consistent with the usual understanding of the deuteron, that the 
tensor force is gained by 3 Si- 3 Di coupling. Thus, we found that, in the description 
of a single nucleon wave function, \<fii), the superposition of more than two Gaussian 
wave packets with independent spin wave functions is effective to gain the tensor 
force in the AMD framework. 

Since the operator <ji ■ r changes the parity of a single-particle state, a parity- 
violated single-particle state is necessary to incorporate the effect of this operator. 
In order to properly take into account the parity violation of single-particle wave 
functions, the parity projection before the variation of the total system is important. 
We show the effect of the parity projection in the deuteron in Table I. In the upper 
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row ("w"), we show the result with the parity projection, which corresponds to 
N = 4 wave packets with different spins in Fig. 2(a). The lower row ("w/o") shows 
the expectation values for the wave function obtained by switching off the parity 
projection. The contribution of the tensor force without the parity projection (—1.5 
MeV) is much smaller than that of the result with the parity projection (—10.1 MeV). 
L 2 is also reduced from 0.52 to 0.05, which means a decrease of the D-state mixing 
by switching off the projection. Thus, we confirm the importance of the parity 
projection. For the deuteron calculated with parity projection, its single nucleon 
wave function is actually a parity- violated state with the mixing of s- and p-orbits, 
as shown in §4. 

3.2. Flexibility on the isospin wave function and charge projection 

We calculate the 4 He by using F-T potential, and investigate the effect of the 
point ii) flexibility of the isospin wave function. The results are summarized in Table 
II. In calculations except for "Usual AMD", each single nucleon wave function is 
described by four wave packets (N = 4) with independent intrinsic spin functions 
(different spins). The top row ("Usual AMD") is obtained by the ordinary version of 
AMD, in which each nucleon is described with a single wave packet, and parity and 
angular-momentum projections are performed. The middle row ("Isospin fixed") 
is obtained by a calculation in which each isospin wave function, \^f a ), is fixed to 
be a proton or neutron, \j a ) = \p) or \n), in the same way as the calculation of the 
deuteron. The bottom row ("Isospin free") is obtained by a calculation in which |7^) 
can describe a proton-neutron mixed state, \j a ) = (§ + 1 a ) \p) + (\ ~ la) l n )> with 
the variational parameter 7^, as shown in Eq. (2-6). By comparing these results, 
\Vt\ and L 2 are larger in the case of "Isospin free" than in the case of "Isospin fixed" . 
Therefore, it is found that the use of the isospin-mixed single nucleon wave function 
is important to gain the tensor force in 4 He. 

In the previous and present subsections, we qualitatively show that two exten- 
sions of the AMD framework, i) a parity- violating mean field and a strong correlation 
between spin and space and ii) flexibility of the isospin wave function, are necessary 
to treat the tensor force. However, the results obtained by the F-T potential still fail 
to quantitatively reproduce the experimental values of such properties as the binding 
energies and the root-mean-square radii of these nuclei. Both of the deuteron and 
4 He are overbound and overshrunk, as can be seen in Table I and II. 

We consider that these failures are caused by the too-large contribution of the 
central term of the F-T potential. Due to the extension of the AMD framework, the 



Table I. Effect of parity projection in the deuteron, which is calculated with the F-T potential. 
The upper (lower) row shows the result obtained with (without) parity projection. T, Vc, Vt 
and B.E. are the kinetic energy, central force, tensor force and total binding energy, respectively. 
They are in units of MeV. L 2 is the expectation value of the squared orbital angular momentum. 



Parity projection 


T 


V c 


V T 


B.E. 


L 2 


w 


24.3 


-18.1 


-10.1 


3.9 


0.52 


w/o 


16.0 


-15.7 


-1.5 


1.2 


0.05 
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contribution of the tensor force becomes large, as we expected. However, as shown 
in "Usual AMD" in Table II, the central force in the F-T potential is so strong that 
4 He is already sufficiently bound without any attraction by the tensor force. As a 
result of the large contribution of the tensor force gained by the extended framework 
in addition to such large contribution of the central force, overbinding and shrinkage 
occur. We can explain these failures as being attributed to a mismatching between 
the extended model space of the present framework and the F-T potential. The 
F-T potential should be used in a simple model space, as written in the original 
paper. 19 ) In a simple model space, an explicit tensor term has only a minor effect, 
while the dominant contribution of the original tensor force is renormalized into the 
large attraction of the central term in the F-T potential. It can well reproduce the 
total binding energy of 4 He when it is used in a small model space. In the case 
of the present calculation, due to the relatively large model space of the extended 
AMD, the contribution of the tensor force becomes large. Consequently, since the 
tensor force is already renormalized into the central force, it is doubly counted in 
the central term and in the direct tensor term, causing such overbinding. 

3.3. Different width wave packets and angular-momentum projection 

As explained in the former subsection, we can not quantitatively reproduce the 
properties of the deuteron and 4 He with the F-T potential, which has a strong 3 E 
central force. Since this potential is not suitable for our research strategy, where we 
treat the tensor force directly without renormalizing into the central force, we employ 
another type of potential based on the AK(T) potential, in which the magnitude of 
the 3 E central force is nearly equal to that of the tensor force. We think that the 
AK(T) potential is more suitable to the extended AMD framework than the F- 
T potential, because the tensor force is not renormalized into the central force. As 
discussed in Ref . [5] , the matrix element of the tensor force is reduced approximately 



by a factor of w ^ in simple mean-field approaches. Therefore, we enhance the 



strength of the tensor force by some factor also in the case of AMD in order to 
properly evaluate the matrix element of the tensor force. In the present calculation, 
we phenomenologically make the strength of the tensor force in the original AK(T) 
potential twice enhanced, so as to reproduce the binding energy and radius of 4 He. 

We calculated 4 He with the AK(T) potential under various conditions of the 
model space and the variation, and investigated the effect of different width wave 



Table II. Results of 4 He with the F-T potential. T, V c , Vls, Vt, V e and B.E. are the kinetic energy, 
central force, LS force, tensor force, Coulomb force and total binding energy, respectively. They 
are in units of MeV. L 2 is the expectation value of the orbital angular momentum. R rm s is the 
root-mean-square radius of the system in units of fm. The experimental values of the binding 
energy and R rm3 of 4 He are B.E. = 28.3 MeV and R rms = 1.47 fm. 23) 




T V c 



V LS 



V e Bl 17 



'rms 



Usual AMD 
Isospin fixed 
Isospin free 



56.0 -85.3 
57.4 -85.7 
77.6 -96.2 



0.0 
0.2 
2.1 



-1.7 
-4.2 
-22.4 



0.9 30.1 0.01 
0.9 31.4 0.06 
1.0 37.9 0.41 



1.37 
1.36 
1.22 
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packets (point iii)) and that of the angular momentum projection (an ingredient 
of the point iv)). We list the conditions I ~ V in Table III. Except for case V, 
the single-nucleon wave functions are represented with four wave packets, which 
corresponds to the case TV" = 4 in Eq. (2-4). We perform the calculation of 4 He with 
a common width parameter for all packets (I, II: z/s="com.") and with different 
width parameters of wave packets (III, IV: z/s="diff."). In the former case "com.", 
a common width parameter v a = v is used, and is chosen to be v = 0.6 so as to 
minimize the energy of 4 He. In the latter case of different widths ("diff."), v a is 
given as 

v a = v x x {vn/vx)^ . (3-1) 

After some trials of various sets {v\, 1/4}, we found the optimum set, v\ = 0.3 and 
= 1.5, which gives the energy minimum. The adopted width parameters {u a } 
are listed in Table III. In each case of a common width (I, II) and different widths 
(III, IV), we first performed the energy variation without the angular-momentum 
(J) projection (I, III), and projected the obtained wave functions onto the J-eigen 
states (II, IV). 

The results of 4 He calculated under each condition are given in Table IV. By 
comparing II (a common width) and IV (different widths), it is found that the 
reproduction of both the binding energy and the radius is improved in result IV. 

Table III. Summary of conditions in calculating 4 He with the AK(T) potential. "JV" means the 
number of wave packets for a single nucleon wave function. V's" indicates that the width 
parameters of the wave packets are equal to ("com.") or different from each other ("diff."). 
Vi ~ v?" show values of the adopted width parameters in units of fm~ 2 . "J pro." indicates 
that the angular momentum projection has been done ("w") or not ("w/o"). "w/o" in "J 
const." means that a solution is obtained by the usual frictional cooling method without a 
J constraint. We constructed a better solution by the frictional cooling method with the J 
constraint ("w" in "J const."). A detailed explanation is given in the text. 





N 


v's 


Vl 




V-6 


v A 


v 5 


J pro. 


J const. 


I 


4 


com. 


0.60 


0.60 


0.60 


0.60 




w/o 


w/o 


II 


4 


com. 


0.60 


0.60 


0.60 


0.60 




w 


w/o 


III 


4 


diff. 


0.30 


0.51 


0.88 


1.50 




w/o 


w/o 


IV 


4 


diff. 


0.30 


0.51 


0.88 


1.50 




w 


w/o 


V 


4 


diff. 


0.30 


0.51 


0.88 


1.50 




w 


w 


V 


5 


diff. 


0.30 


0.45 


0.67 


1.00 


1.50 


w 


w 



Table IV. Results of 4 He with the AK(T) potential calculated under various conditions. The con- 
ditions "I ~ V " correspond to TABLE III. J 2 , L 2 and S 2 are the expectation values of the 
total angular momentum, the orbital angular momentum and the spin, respectively. 





T 


V c 


V T 




B.E. 


Rrms 


J 2 


L 2 


S 2 


I 


94.1 


-39.4 


-66.4 


0.9 


10.8 


1.23 


0.32 


1.50 


1.31 


II 


91.7 


-40.3 


-68.0 


0.9 


15.7 


1.23 


0.00 


1.27 


1.27 


III 


79.0 


-41.9 


-55.4 


0.9 


17.4 


1.32 


0.23 


0.80 


0.68 


IV 


73.4 


-43.5 


-56.2 


0.9 


25.3 


1.32 


0.00 


0.57 


0.57 


V 


70.4 


-43.5 


-56.3 


0.9 


28.6 


1.35 


0.00 


0.50 


0.51 


V 


69.0 


-43.0 


-55.5 


0.9 


28.7 


1.35 


0.00 


0.54 


0.54 
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This shows the importance of using wave packets with different widths. 

The comparison of III (without J projection) and IV (with J projection) indi- 
cates the importance of the angular-momentum projection. As can be seen in Table 
IV, the binding energy gains 8 MeV in result IV by the J projection. Before the 
J projection, the J 2 value is small, as can be seen in result III, which indicates 
that the J / components are very small. Considering that the state is almost the 
J = eigen state already before the J projection, the energy gain by the projection 
seems to be unexpectedly large. This gain is dominantly attributed to a decrease of 
the kinetic energy term. Dissolving the III state (|<£m)) into various (J, J z ) compo- 
nents (| J J z )), we analyzed the contributions of each (J, J z ) component to the kinetic 
energy (Table V). Detailed expressions of |<Pm) and \ J J z ) are 

l^m) = P Tz =oPp=+m, (3-2) 
\JJz) = Pj,j z \$m), (3-3) 

where |^) is the intrinsic state of the III state. As shown in Table V, the J = 
component is dominant as 96%, while the J = 1 and J = 2 components are mixed 
into the III state at the rate of a few percent. Although the J = 1 and 2 components 
are very small, they have huge kinetic energy (~ 200 MeV), and hence the significant 
gain of the kinetic energy is caused by excluding the J / components after the 
J projection. Such a huge kinetic energy of the J / components arise from very 
narrow wave packets, such as v± = 1.5 used in calculating III and IV. It is considered 
that high-momentum components are contained in the J ^ components of the III 
state, due to such narrow wave packets. 

Since our calculation procedure is "Variation Before Projection" with respect to 
the J projection, we can expect to find a better solution of the J = state than IV. 
For this aim, we apply a constraint condition method in the present framework. We 

* 2 

put a constraint on the squared total angular momentum by setting W2 = ( J ) in the 
Eq. (2-11) for the energy variation, and constructed the states with various values 
of (J ) (J is the total angular momentum operator). After the energy variation, we 

Table V. Contributions from various J components to the kinetic energy. "Before J projection" 
corresponds to the III state in Table IV. "overlap" indicates the ratio of each (J, J z ) component 
(| J J z )) included in the III state (|*in>); | ( J Jz |#m) j 2 /VV Jz I J Jz) (^iiil^m). The definitions 
of I $111} and I J J z ) are explained in the text. 



J 


Jz 


T 


overlap 


T x overlap 


Before 


J projection 


79.04 








73.43 


0.955 


70.09 


1 


1 


199.75 


0.007 


1.38 


1 







0.000 


0.00 


1 


-1 


201.24 


0.007 


1.35 


2 


2 


238.19 


0.005 


1.23 


2 


1 


172.96 


0.002 


0.31 


2 







0.000 


0.00 


2 


-1 


177.75 


0.002 


0.30 


2 


-2 


238.38 


0.005 


1.24 
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projected them onto the J = state and obtained an energy surface as a function 

* 2 

of the value of the constraint (J ). We then chose the energetically best solution, 
which is state "V" shown in Table IV. In this solution, the binding energy and the 
radius are more improved than IV. In this final result, V, the binding energy and 
the radius are —28.6 MeV and 1.35 fm, respectively. We checked the convergence 
of the solution with an increase of the number of wave packets, N. V is the result 
obtained by using five wave packets for each nucleon wave function, namely N = 5 
in Eq. (2-4), in the same way as the case of V for N = 4. Since the changes of 
various properties from the V state to the V state are very small, we consider that 
the convergence of the solution is sufficient at N = 4. 

Thus, we confirmed that extension iii) different width wave packets and especially 
the angular-momentum projection in iv) are effective as well as extensions i) and ii). 

3.4. Systematics of deuteron, triton and alpha with AK(T) potential 

Within the extended AMD framework we could rather quantitatively reproduce 
the properties of 4 He by using the AK(T) potential with a twice-enhanced strength 
of the tensor term. We now consider our calculation of the deuteron and triton in 
the same way as in case V for 4 He. The results are given in Table VI. Comparing 
the theoretical values of binding energies and radii with the experimental data, we 
can tolerate the result of the triton, though it is somewhat overbound, whereas the 
deuteron is terribly overbound and overshrunk. 

One of the reasons for this failure in systematic reproduction is conjectured to 
be as follows: The deuteron is a special system, because it is composed of only two 
nucleons. It has just one degree of freedom in the spatial coordinate, namely, the 
relative motion between them. Since its relative motion may be described sufficiently 
well by the superposition of Gaussian wave packets in the extended AMD, the two- 
body correlation in the deuteron is incorporated much better than that in 4 He. As 
mentioned in the previous section, Ref. [5] pointed out that the matrix element of 

the tensor force is reduced approximately by a factor of ^J~^ in the simple mean-field 

type treatment of 4 He due to the limitation of the model space, which is insufficient 
to incorporate a two-body correlation. In order to regain this reduction, we need a 



Table VI. Results of the deuteron, triton and 4 He calculated with AK(T) potential. "N" means 
the number of wave packets for a single nucleon wave function, "it" indicates the enhancement 
factor of the strength of the tensor force, "ai", "b" and "c" show the results of the deuteron, 
triton and 4 He calculated with the twice-enhanced tensor force. "a2" is the deuteron calculated 
with a 1.5-times enhanced tensor force. The experimental values of the binding energy and 
radius: B.E. = 2.2 MeV and R rms = 1-94 fm for deuteron, and 8.5 MeV and R rms = 1-55 fm 
for triton. 23) 





N 


X T 


T 


V c 


V T 


v e 


B.E. 


Rrms 


J 2 


L 1 


S 1 


ai 


2 R 


5 


2.0 


34.9 


-4.1 


-44.1 


0.0 


13.4 


1.00 


2.00 


0.72 


1.98 


a2 


2 H 


5 


1.5 


23.2 


-4.2 


-20.7 


0.0 


1.7 


1.19 


2.00 


0.49 


1.99 


b 


3 H 


5 


2.0 


48.0 


-20.5 


-38.0 


0.0 


10.5 


1.27 


0.75 


0.54 


1.01 


c 


4 He 


4 


2.0 


70.4 


-43.5 


-56.3 


0.9 


28.6 


1.35 


0.00 


0.50 


0.51 
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twice-stronger tensor force for the calculation of 4 He. However, because of the better 
description of the two-body correlation in the deuteron wave function than in 4 He, 
the reduction factor of the matrix element must be smaller and the enhancement 
factor of the strength should be smaller in the deuteron than in 4 He. This means 
that the twice-enhanced tensor force is too strong for the deuteron. Actually, from 
the value L 2 = 0.72 shown in row "ai", we can find that the -D-state probability of 
the deuteron is too large (about 12%), compared with the well-known value of several 
%. We then calculated the deuteron with a smaller enhancement factor, xt = 1-5, 
of the strength of the tensor force. The result is shown in row "a2" of Table VI. In 
this case, L 2 = 0.49, which indicates that the D-state probability is estimated to be 
a reasonable value as 8%. The binding energy also decreases to be 1.7 MeV, which 
well agrees with the experimental data. These results support that an about 1.5 
time enhancement of the tensor force is suitable to reproduce structure properties in 
the case of the deuteron system. 

§4. Analysis of and discussion on single particle levels 

The 4 He obtained by the extended version of AMD with the AK(T) potential 
has a significant contribution of the tensor force to the potential energy. Due to the 
tensor force, the obtained wave function is expected to be more complicated than 
the ordinary (0s) 4 -like wave function obtained with the usual effective forces with 
no tensor term. In this section, we clarify its internal structure. 

In the procedure of the present calculation to obtain a total wave function, 
we operate parity-, charge- and angular- momentum-projections to an intrinsic wave 
function. The wave function of 4 He is written as 

| 4 He) = Pj = oP Tz =oPr=+\$), (4-1) 

where the intrinsic state, \<P), is expressed by a single Slater determinant. We can 
analyze the Hartree-Fock-like single-particle orbits, \ai) (i = 1 ~ 4), (AMD-HF 
orbits) extracted from the AMD wave function, |#), by the method given in Ref. 
[22]- ' 

Table VII shows the properties of the AMD-HF levels extracted from the intrinsic 
wave function of result V in Table IV. Apparently, the intrinsic state \<P) is not the 
simple (0s) 4 state. If |#) were the (0s) 4 state, j , I 2 and P(+) should be equal 

Table VII. Single particle levels of 4 He in the state V. "S.P.E." is the single particle energy of each 
level in unit of MeV. "j 2 " and "i 2 " are the expectation value of total angular momentum and 
orbital angular momentum, respectively. "P(+) (P(— ))" and "P(proton) (P(neutron))" are the 
ratio of the positive(negative) parity component and the proton(neutron) component included 
in each level, respectively. They are in unit of %. 



Level 


S.P.E. 


3 2 




P(+) 


PR 


P (proton) 


P(neutron) 


1 


-4.04 


1.16 


0.83 


81 


19 


61 


39 


2 


-4.52 


1.17 


0.85 


81 


19 


59 


41 


3 


-8.28 


1.09 


0.68 


85 


15 


30 


70 


4 


-8.78 


1.06 


0.65 


86 


14 


45 


55 
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to 0.25, 0.00 and 100% at each level, respectively. The obtained values, however, 
deviate from these values. The obtained levels contain a negative-parity component 
with 15 ~ 20%, and I 2 is non-zero. These results indicate a mixture of the p state 
in all single-particle orbits. It is interesting that they are also proton-neutron-mixed 
states, as shown in Table VII. 

To investigate in more detail, we have decomposed the AMD-HF level into the 
eigen-state of parity and charge. The composition of each level is shown in Table 
VIII. I 2 and the root -mean-square radius of each component are given in Table IX 
and X, respectively. The negative-parity components included in all AMD-HF levels 
are found to be dominated by p orbits, because I 2 for the negative-parity component 
for both of proton and neutron is nearly equal to 2, as shown in Table IX. In the 
positive-parity components, non-zero I 2 values imply the mixing of d orbits with the 
dominant s orbit. As shown in Table X, the radius of the negative-parity components 
is about 1.28 fm. If we assume a harmonic-oscillator wave function, this radius, 1.28 
fm, of the p-orbits suggests an unusually small oscillator size. On the other hand, the 
radius of the positive-parity components is about 1.6 fm, which suggests a normal 
size for the s obit in 4 He. Thus, it is considered that the narrow p state (negative- 
parity component) is mixed with the normal-size s state (positive-parity component) 
in each level. 

In order to confirm the above picture, we analyze the AMD-HF orbits with 



Table VIII. Composition of the single-particle levels, proton (neutron) + ( — ) means the component 
of the proton (neutron) with the positive (negative) parity. The values are in unit of %. 

Level proton + neutron + proton — neutron — 



1 


55 


26 


6 


13 


2 


50 


30 


9 


10 


3 


21 


64 


9 


6 


4 


37 


49 


8 


6 



Table IX. Expectation value of I 2 of each component. 



Level proton + neutron + proton — neutron — 



1 


0.16 


0.99 


2.58 


2.56 


2 


0.29 


0.64 


2.61 


2.61 


3 


1.09 


0.10 


2.68 


2.63 


4 


0.41 


0.26 


2.67 


2.48 



Table X. Root-mean-square radius of each component. The units are fm. 



Level proton + neutron + proton — neutron — 



1 


1.59 


1.70 


1.30 


1.27 


2 


1.58 


1.60 


1.27 


1.31 


3 


1.55 


1.54 


1.26 


1.28 


4 


1.52 


1.51 


1.28 


1.25 
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Fig. 3. Overlaps between the AMD-HF states extracted from state V in Table IV and the s and 
p states in the harmonic-oscillator potential with various width parameters {v s and v p fm~ 2 ). 
(a) Overlap between the positive parity component of each AMD-HF state and the s state; 
|(sH.o.|a^)j 2 - (b) Overlap between the negative parity component of each AMD-HF state and 
the p state; |(pH.o.|a~)| 2 . The four lines in each figure correspond to four single-particle levels, 
respectively. 




0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 



u s [fm" 2 ] u p [fm- 2 ] 

(a) positive parity (b) negative parity 



simple s- and p-orbits written in terms of harmonic-oscillator (H.O.) wave functions: 



1 



rn 



\sh.o.) ■ (r\0si,m,T z )= R s (r 

R s (r) = l\:~ 
Iph.o.) : (r\0pj,m, t z ) = R p (r 



1 



exp [-v s r 2 ] , 
\^m h -m s \jm\Y lmi {Q) 



.mi,m s 
l 



,m.< 



R p (r) = N p 2 r exp [-u p 



(4-2) 
(4-3) 

(4-5) 



Here, v s and u p indicate the width parameters of the s and p orbits, respectively. 
Fig. 3(a) (3(b)) shows the overlap between the positive (negative)-parity component 
of the AMD-HF states, \af), and the sh.o. (Ph.o.) state as a function of the H.O. 
width parameter, v s (f p ). Here, \af) is the parity-eigen state projected from \cti), 



Table XI. Composition of the p state at v p = 0.8 fm -2 . The units are %. "p" and "n" in the 
column "p or n" mean the proton and neutron component, respectively, "p+n" means the sum 
of proton and neutron components. 



p or n 


level 1 level 2 level 3 level 4 


Opi P 
Opi n 

2 

Opi p+n 


27.2 41.3 54.8 48.0 
59.8 43.6 32.4 39.2 
87.0 84.9 87.2 87.3 


Op 3 p 

2 

Op 3 n 

2 

Op 3 p+n 

2 


2.0 2.9 3.3 3.6 
3.4 4.0 1.8 2.2 
5.4 6.9 5.1 5.8 


0p p+n 


92.4 91.8 92.3 93.1 
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Fig. 4. Comparison of the radial dependence of the single-particle wave functions and the tensor 
force. "R s * R p " (dashed line) is the product of the radial wave function of the s state (Eq. 4-3) 
and that of the p state (Eq. 4-5). "Vr" (solid line) is the radial part of the tensor force of the 
AK(T) potential. Both are depicted after being multiplied by r 2 . 




_g I 1 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 4 
r[fm] 



where {af\af) = 1. |sh.o.) and \pn.o.) have maximum overlap (92%) with \af) 
and \a~) at v s = 0.3 fm~ 2 and u p = 0.8 fm -2 , respectively. On the other hand, 
the typical width parameter, v = 0.26 fm~ 2 , is known to reproduce the size of 4 He 
by the H.O. (0s) 4 wave function. Thus, it is explicitly confirmed that the narrow 
p state is mixed with the normal-size s state. It is worth showing that the 92% Op 
state consists of 86% Opi and 6% 0^3 (Table XI). 

2 2 

The reason for the mixture of such a shrunk p orbit can be understood as follows. 
In the single-particle picture, the matrix element of the tensor force is gained by 
the coupling of the s state and the p state, such as ( [pp]d \vj 2 \ [ss]s ) in which two 
nucleons in the p orbit and those in the s orbit are coupled to the d orbit and the 
s orbit, respectively. Fig. 4 shows the product of the radial wave functions of the 
•sh.o. orbit (y s = 0.3) and the narrow pn.o. orbit {y v = 0.8), and also the radial 
part of the tensor force in the AK(T) potential. The peak position of the product 
of these two orbits is consistent with the potential range of the tensor force. Such 
a combination as sh. 0.(^3 = 0.3) and ph.oX^p = 0.8) is suitable to incorporate the 
tensor effect into the binding energy. As a result of the analysis of the single particle 
levels, we are convinced of the role of extension iii). In other words, we should use 
wave packets with different width parameters in order to describe the shrunk p orbit 
simultaneously with the normal-size s orbit. 

§5. Summary 

In the present paper, we tried to construct a new model that can treat the tensor 
force directly, based on the antisymmetrized molecular dynamics method (AMD). 
We employed the Furutani-Tamagaki potential (F-T potential) and the Akaishi(T) 
potential (AK(T) potential). As a result of practical calculations, we confirmed 
that four extensions in the AMD framework are effective for a direct treatment of 
the tensor force: i) a parity-violating mean field and a strong correlation between 
spin and space, ii) flexibility on the isospin wave function, iii) different width wave 
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packets and iv) charge-, parity- and angular-momentum-projections. Due to all of 
the above extensions, we obtained a satisfactory solution of 4 He. We can reproduce 
the binding energy for 4 He with the AK(T) potential with a twice-enhanced tensor 
strength; also, the radius of the calculated 4 He well agrees with the experimental 
value; The binding energy and root-mean-square radius were calculated to be 28.6 
MeV and 1.35 fm, respectively. The energy composition is as follows: T = 70.4 
MeV, V c = -43.5 MeV, V T = -56.3 MeV and V Co ulomb = 0.9 MeV. Thus, this 
solution has a large contribution of the tensor force, as we expected. We found that 
the influence of the large contribution of the tensor force appeared in single-particle 
wave functions of this 4 He solution. The shrunk p orbit was mixed to the usual 0s 
orbit to gain the tensor force. 

We could gain a large contribution of tensor force in the deuteron, triton and 
4 He with the extended version of AMD. However, we are not satisfied with the 
results because we need a two- fold enhancement of the tensor force in the AK(T) 
potential to quantitatively describe the properties of 4 He, and have not yet succeeded 
to reproduce the systematics of the binding energies of the deuteron, triton, and 
4 He. These fact may mean that the extension of the framework is still insufficient. 
Referring to the deuteron system, we here comment on possible residual correlations. 
Since the deuteron has only one degree of freedom, i.e. the relative motion, the 
description of the two-body correlation is considered to be directly improved by 
increasing the number of wave-packets. In fact, we have found that the tensor 
contribution in the deuteron changes drastically with an increase of the number 
of wave packets from N = 5 to 10, with its binding energy changing from 13.4 
MeV to 18.2 MeV. It is different from the case of 4 He, where the results almost 
converge at N = 4. In the case of N = 10 for the deuteron, only a 1.4-times 
enhancement of the tensor force is sufficient to get the satisfactory result of 1.6 
MeV for the deuteron binding energy, which agrees with the experimental data (2.2 
MeV). This result implies that there is some scope to improve the description of 
the two-body correlation in our framework. In a mean-field- type model, such as the 
AMD method, it is difficult to treat the short-range part of the tensor force. This 
difficulty is expected to be overcome by introducing the so-called T-type base, 24 ) 
or renormalizing such a part of the tensor force into the central force. In order to 
perform further studies based on AMD, the latter way is suitable for our method. 
Due to renormalization of the short-range tensor force, the central force is expected 
to become more attractive than the present one. From the result obtained with the 
F-T potential, which has a strong central force, we can conjecture that the binding 
energy is easily reproduced by the renormalized force without the enhancement factor 
of the tensor force. In a study with a mean-field approach by Sugimoto et a/., 12 ) 
they use the enhancement factor 1.5 for the same tensor force as in the present 
study. By adjusting the strength of the central force, the binding energy of 4 He 
is well described in their study. Although the central force is phenomenologically 
fitted, their approach is interpreted to simulate the above-mentioned renormalization 
of the short-range tensor force in the central force. Needless to say, it is important to 
derive effective nuclear interactions with a direct long-range part and a renormalized 
short-range part of the tensor force based on realistic forces, as much as possible. 
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Such a renormalization framework may work effectively to investigate the influence 
of the tensor force on the nuclear structure, because the short-range part of tensor 
force is expected to affect mainly the bulk properties of nuclei, while its long-range 
part should be sensitive to the details of the nuclear structure. 
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Appendix 

Parameters in the Akaishi(T) potential are given in Table XII. C„ E , Cf t E , 
and Cn° are parameters of the central force in Eq. (2-2) and T% E is that of the 
tensor force in Eq. (2-3). Each range parameter, b n , is generated from b\ and &io by 
b n = b 1 x(b w /b 1 )( n - 1 )/v. 



Table XII. Range parameters {b n } and coefficients {C*} in vfj and {T„} in vjj of AK(T) potential. 
The units of {&„}, {C*} and {T„} are fm, MeV and MeV/fm 2 , respectively. The definitions 
are given in Eqs. (2-2) and (2-3). 



n 


b„ 










1 n 


1 


0.1600 


-3.4970 xlO 2 


-3.2879 xlO 2 


2.7042 xlO 1 


1.5355xl0 2 


-2.2053xl0 2 


2 


0.2141 


1.4440 xlO 3 


1.2518xl0 3 


-2.0400 xlO 2 


-1.2405 xlO 2 


1.1121xl0 3 


3 


0.2865 


-2.9647xl0 3 


-2.5066 xlO 3 


9.0073 xlO 2 


-5.1523xl0 2 


-1.6079xl0 3 


4 


0.3833 


1.4499xl0 3 


1.3850xl0 3 


-1.8442xl0 3 


4.7590xl0 2 


1.3446xl0 2 


5 


0.5129 


1.6601 xlO 3 


1.2147xl0 3 


1.3668 xlO 3 


4.3360xl0 2 


-3.8447xl0 2 


6 


0.6863 


-7.7942 xlO 2 


-6.5449 xlO 2 


1.8709 xlO 1 


2.5416X10 1 


-1.8164xl0 2 


7 


0.9183 


-1.4783X10 1 


7.6192X10 1 


8.5685X10 1 


-4.6085 xlO 1 


-4.1188X10 1 


8 


1.229 


-6.1748X10 1 


-6.5063X10 1 


-3.5048X10" 1 


4.3151 


-1.2277X10 1 


9 


1.644 


6.5369 


1.1420X10 1 


-6.3206 


-6.8892 


-2.2104 


10 


2.200 


-6.4642 


-6.6031 


1.3007X10 1 


1.6804 


-6.5291 xlO" 1 
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